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Summary. Forty years after the discovery of rotation-powered pulsars, we still do 
not understand many aspects of their pulsed emission. In the last few years there 
have been some fundamental developments in acceleration and emission models. I 
will review both the basic physics of the models as well as the latest developments 
in understanding the high-energy emission of rotation-powered pulsars, with partic- 
ular emphasis on the polar-cap and slot-gap models. Special and general relativistic 
effects play important roles in pulsar emission, from inertial frame-dragging near 
the stellar surface to aberration, time-of-flight and retardation of the magnetic field 
near the light cylinder. Understanding how these effects determine what we observe 
at different wavelengths is critical to unraveling the emission physics. I will discuss 
how current and future X-ray and gamma-ray detectors can test the predictions of 
these models. 



1 Introduction 

Rotation-powered pulsars are fascinating astrophysical sources and excellent 
laboratories for study of fundamental physics of strong gravity, strong mag- 
netic fields, high densities and relativity. The major advantage we have in 
studying pulsars is that we know they are rotating neutron stars and that 
they derive their power from rotational energy loss. The challenge is then to 
understand how they convert this source of power into the visible radiation. It 
is generally agreed that this occurs through acceleration of charged particles 
to extremely relativistic energies, using the rotating magnetic field as a unipo- 
lar inductor to create very high electric potentials. Beyond this fundamental, 
there is a large divergence of thought on what comes next: whether the ac- 
celeration occurs in the strong field near the neutron star surface or in the 
outer magnetosphere near the speed of light cylinder, or even beyond the light 
cylinder in the wind zone. The particle acceleration may well be occurring in 
all of these regions, either in the same pulsar or in pulsars of different ages. 

In recent years, there has been much activity both in new detections and in 
theoretical study of rotation-powered pulsars. Multibeam radio surveys at the 
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Parkes Telescope [S31[5S] have increased the population of known radio pulsars 
to more than 1700. In addition, extended radio observations of supernova 
remnants and unidentified 7-ray sources have discovered a number of young 
pulsars that are too radio- faint to be detected by surveys [13] . Although pulsed 
emission at other wavelengths has been detected from only a small fraction of 
these, this number is growing as well. At the present time, there are 7 pulsars 
with high-confidence detection of 7-ray pulsations [44] , about 30 having X-ray 
pulsations [45] and 10 with optical pulsations [55] . 

This paper will review both the fundamental physics and latest theoretical 
developments of acceleration and radiation in polar cap models. A comple- 
mentary paper by Cheng [14] reviews acceleration and radiation in the outer 
gap model. 



2 Acceleration Near the Polar Cap and Beyond 

A magnetic dipole rotating in vacuum will induce an electric field both along 
and across the magnetic field lines [2^. In the case of a pulsar, with high 
angular velocity Q and surface dipole fields Bq ~ 10^^ G, the electric force 
parallel to the magnetic field just above the neutron star surface exceeds the 
gravitational force by many orders of magnitude. Vacuum conditions therefore 
cannot exist outside a pulsar, since charges can be pulled from the stellar 
surface [26]. If the charge density reaches the Goldreich- Julian value, 

V-E n-B 



derived from the condition 



(O X r) X B 

E = -^ , (2) 



then the electric field parallel to the magnetic field vanishes. This is the force- 
free solution where charges and magnetic field corotate with the star. Coro- 
tation must break down near the light cylinder, Rl = fi/c, due to particle 
inertia. But if vacuum cannot surround a pulsar, neither can a completely 
force-free magnetosphere, since in that case no acceleration of charge, currents 
or radiation would exist. A real pulsar must operate somewhere between the 
two extremes of the vacuum and the force-free states, but a self-consistent, 
global solution has not yet been found. Global magnetospheric simulations 
study how a rotating neutron star magnetosphere fills with charge from the 
vacuum state. The resulting domes and torii of charge that build up near the 
pole and equator [JT] [SS] seem to be unstable [70], but the particle-in-cell 
codes cannot run long enough to reach a stable magnetospheric configuration. 
However, it seems that a near force-free magnetosphere |71| cannot result only 
from charges flowing out of the stellar surface, but requires some extra source 
of charge created in the magnetosphere above the surface. This extra source of 
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Fig. 1. Illustration of space-charge limited flow and vacuum gap accelerators above 
a pulsar polar cap. Ts is the neutron star surface temperature and Ti,e are the ion 
or electron thermionic temperatures. 



charge is thought to be production of electron-positron pairs by the photons 
radiated by accelerating particles. The pulsar magnetosphere must somehow 
be made up of self-consistent force-free and non-force-free regions in balance 
with each other. One way to determine the structure of these regions is to 
study the microphysics of the electrodynamics and charge flow at different 
sites where acceleration may occur. One of these sites, the polar cap accelera- 
tor, is near the neutron star surface on the open field lines that cross the light 
cylinder. 



2.1 Polar cap accelerators 

The two main types of polar cap accelerator are vacuum gaps [671 176] and 
space-charge limited flow (SCLF) gaps [6l [35]. There are binding (or cohe- 
sive) forces on charged particles in the neutron star surface due to the lattice 
structure in a strong magnetic field, such that particles are free only if the 
surface layers are above the thermionic emission temperature. For electrons, 
this temperature is [75] 

T..3,.xlO»A-(|)°"(^)" (3) 



and for ions it is 
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R \ 0.73 



T.^3.5xl0^i.^^J (4) 

where Bq is the surface magnetic field strength and Z is the atomic num- 
ber of matter in the surface layer. If the surface temperature, < Ti^e, 
then charges are trapped in the surface and the full vacuum electric field, 
E\\ = i7i?o^, where R is the neutron star radius, exists above the surface. 
If Ts > Ti^e, then charges are "boiled off' the surface layers and can flow 
along the open field lines in SCLF. Measured surface temperatures of pulsars 
are typically Tg > 0.5 — 3.0 x 10^ K, above Te and Ti for the normal range 
Bq < 10^'^ of surface fields. However, a few high-field pulsars or magnetars may 
have Ts < Ti^e and vacuum gaps [76l [81] . If > Ti^e and the full Goldreich- 
Julian charge can be supplied right about the surface, the true charge density 
along each field line will drop faster, p oc r~^, than poj as the dipole field 
lines fiare. A SCLF electric field, V • ~ {p — Pgj)/^o grows with altitude. 
The two types of accelerator thus differ by the surface boundary condition: 
where p{R) = 0, ^ for vacuum gaps and p{R) = pcj, E\\{R) = 

for SCLF accelerators. The form of in SCLF accelerators is thus sensitive 
to the detailed distribution of the charge density, which depends both on the 
open field line geometry as well as the compactness of the neutron star. At 
altitudes {z <C 9pc, z :s> dpc}, with z = {r/R — 1) being the height above 
the surface. 



{^'^Pc(i) ^}'^cosa + {z,el^ (j^) ^^^}%^sinacos^ 



-E'li — Bq6 

|621 135] . where 9 and ip are the magnetic polar and azimuth angles, a is the 
magnetic inclination angle to the rotation axis, k = 2GI/{c^R^) is the stellar 
compactness parameter, 9^^ ~ {QR/cY^^ is the polar cap half-angle and / 
the neutron star moment of inertia. The first term in Eqn ([5]) is due to inertial 
frame dragging near the neutron star surface, and dominates for small r and 
low inclination, while the second term is due to the flaring of the field lines. 

Both accelerators will be self-limited by the development of pair cascades, 
at altitudes where the particles reach high enough Lorentz factors to radi- 
ate 7-ray photons. The dominant process in most pulsars is one-photon pair 
production, which can occur only in very strong magnetic (or electric) fields. 
In this process, the magnetic field absorbs the extra momentum of a photon 
having the energy required to create a pair, so that the threshold condition on 
the photon energy is = 2mc^ / sin9jB, where 9jb is the angle between the 
photon propagation direction and the local magnetic field. The accelerated 
particles moving along magnetic field lines with high Lorentz factors radiate 
7-ray photons at very small angles to the field {9o ~ I/7) , so the one-photon 
pair production rate for these photons is initially zero. However, as they prop- 
agate through the curved dipole field, their angle increases until the threshold 
condition and the attenuation coefficient becomes large. The voltage across a 
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vacuum gap breaks down when a stray 7-ray crosses the magnetic field within 
the gap and creates a pair. The electron accelerates upward/downward for 
fl ■ B{>, <)0 over the polar cap, and the positron accelerates in the opposite 
direction as shown in Figure [1] Both particles produce more pairs when their 
radiated photons reach the pair threshold, causing a pair avalanche and sud- 
den discharge of the vacuum. The potential drop in the gap thus oscillates 
between Kg ~ f2Bo[Re^^f/2 and 0. In SCLF for [2 ■ B(>,<)0, an elec- 
tron/positron accelerates upward from the surface until the radiated photons 
reach pair threshold, at which point the positrons/electrons from the pair 
decelerate, turn around, and accelerate downward toward the neutron star 
surface. The polarization of pairs above the pair formation front (PFF), dis- 
cussed in more detail in Section [3l may short out the E\\, halting any further 
acceleration at higher altitude. These accelerators can thus maintain a steady 
current of upwardly accelerating electrons, at ~ cRgj, and a downward 

current of positrons, at <JC cpcj, which heat the polar cap. The accelerator 
voltage is determined by the height of the PFF, which is again roughly com- 
parable to the pair creation mean-free path. However, the stability of SCLF 
accelerators has not yet been verified through time-dependent models, and 
some simplified studies [49] have in fact shown that some oscillations in the 
pair creation rate could exist. This is certainly an issue that needs further 
investigation. Stability of vacuum gaps has also been studied by Gil et al. 

m- 

2.2 Death lines 

The pair cascades can be initiated either by curvature radiation (CR) |16j 
or by resonant or non-resonant inverse-Compton scattering (ICS) of stellar 
thermal X-rays by primary electrons [72| . Since for a given Lorentz factor the 
peak CR photon energy, — 3Xcj^ /2pc, where Ac = h/mc is the electron 
Compton wavelength and pc is the magnetic field radius of curvature, is much 
lower than the ICS peak energy, e-y ~ 7 in the extreme Klein-Nishina limit, 
pair production of CR photons requires a much higher particle Lorentz factor. 
The PFF for ICS therefore occurs at a lower altitude than the PFF for CR 
[35] , The PFF altitude is the sum of the acceleration length oc and the 

pair attenuation length, which arc both inverse functions of the magnetic field, 
and effectively of the pulsar age. If the PFF is larger than a stellar radius, 
then the magnetic field becomes too weak for pair creation to occur and a 
PFF does not exist. For SCLF accelerators, CR photons can produce pairs 
only in the case of young pulsars (r < 10^ yr) and a few millisecond pulsars 
[36l|40j. One can define a 'death line' in P — P space, where P is the period 
derivative, shown in Figure [2] below which pulsars cannot produce pairs from 
CR photons. Below the CR pair death line, pulsars can produce pairs only 
from ICS photons [37]. Below a lower, ICS pair death line, pulsars cannot 
produce any pairs and are expected to be radio quiet. In fact the predicted 
ICS death line falls close to the edge of the known radio pulsar population, as 
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Fig. 2. Period derivative vs. period for radio and high-energy pulsars from the 
ATNF catalog, showing death lines [Sj for production of electron-positron pairs 
through curvature radiation (CR) and inverse Compton scattering (ICS). The ICS 
pair death line was determined for a neutron star surface temperature of lO" K and 
a standard equation of state [8]. 

shown in Figure[2l[38j. The few millisecond pulsars that lie below the ICS pair 
death line may be able to produce pairs through interaction of ICS photons 
with thermal X-ray photons from a hot polar cap [83] . 

3 Electric Field Screening and Polar Cap Heating 

In SCLF accelerators, the polarization of charge above the PFF acts both 
to screen the and to produce heating of the polar cap by the downward 
flowing particles [4J . Figure [3] illustrates the dynamics of electric field screen- 
ing. Primary electrons (e~) accelerate upward from the stellar surface and 
produce pairs at different altitudes above the PFF. The positrons decelerate 
and turn around in a distance short compared to the PFF altitude and each 
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Fig. 3. Illustration of electric field screening in SCLF models. e~ represents the 
primary electron flux accelerating upward from the neutron star surface, e"*" and 
e~ represent electron-positron pairs created by the primary electron above the pair 
formation front. 



reversing positron creates a small excess of negative charge. As more positrons 
are produced and decelerated, the space charge becomes more negative until 
the entire charge deficit 5p = {p — pgj) that produced the £'|| is accounted for. 
Since the charge deficit is small compared to the primary charge (dp Pgj), 
the screening length scale is a very small fraction of the PFF altitude. The 
flux of returning positrons, as a fraction of the primary flux, is approximately 



PGJ 



'2pGJ 



3 

20" 



(6) 



where is the PFF height above the stellar surface and k, is the general 
relativistic factor appearing in the SCLF electric field (Eqn(l5])). The corre- 
sponding polar cap heating luminosity is 



T max 



(7) 



where ' 



is the primary particle flux and ^{zq) is the potential drop at zq. 



The growth in the charge density above the CR PFF is very rapid since the CR 
peak energy e oc 7'^ increases rapidly with increasing electron Lorentz factor, 
7. The CR initiated cascades have very high multiplicities and the screening 
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of ii'ii takes place in a relatively short distance above the PFF [35]. The ICS 
initiated cascades on the other hand have much lower multiplicities since the 
ICS photon production rate decreases with increasing 7, as ~ I/7, for non- 
resonant scattering and first increases sharply, then decreases as ~ 1/7^ for 
resonant scattering. The screening above ICS PFFs takes place over a larger 
scale length, depending on the polar cap temperature and for a number of pul- 
sars well below the CR pair death line, there is no screening [S7\ . Even when 
ICS screening is locally complete, an unscreened charge deficit can develop 
at higher altitudes since the charge density deficit 5p grows faster than the 
increase in pair density. In that case, an reappears and acceleration con- 
tinues. If a CR PFF forms at higher altitude, then screening occurs, otherwise 
unscreened acceleration continues to high altitude [S^ . 

The CR pair heating luminosity is much higher than the ICS heating lu- 
minosity because the CR PFF occurs at higher altitude and a larger flux of 
positrons returns through a higher voltage, bombarding the polar cap with 
higher energy. Figure [4] shows calculated heating luminosities from CR and 
ICS positrons as a function of pulsar characteristic age, Tgd = -P/2P, for sev- 
eral different periods and a PC surface temperature of T = 3 x 10^ K. The 
CR heating lines terminate on the high Tgd side at the age corresponding to 
the CR pair death fine, which is around Tgd ~ 10^ yr for periods P = 0.1,0.2 
s and around Tgd ~ 10* yr for periods P = 5 ms. The CR heating luminosities 
are several orders of magnitude higher than the ICS heating luminosities, for 
normal pulsars. All the L+ are higher for millisecond periods, because the 
PFFs occur at higher altitude and the gap voltage is larger. A sudden drop in 
L+ is therefore predicted at the CR death line. Pulsed X-ray emission has now 
been detected from many middle-aged and older pulsars [9]. Several thermal 
components are often seen, one having lower T and larger area, which may be 
full surface cooling, and another having a higher T and smaller area, which 
could be due to polar cap heating. The luminosity of the hot thermal compo- 
nents, which are shown in Figure [H may indicate that the heating luminosity 
drops suddenly around Tsd ^ 10'' yr or somewhat sooner. To within an order 
of magnitude, the observed luminosities of the hot thermal components agree 
with the predicted L+, for both normal and millisecond pulsars. 

Although the i+ plotted in Figure [4] were computed numerically, analytic 
expressions for L+ provide a good estimate in the case of complete screening 
and upper limits in the case of incomplete screening. The heating from CR 
pair fronts predicts a surface X-ray luminosity of approximately [36] 



where B0.12 = Bq/IO^^ G. The heating from ICS pair fronts predicts a surface 
X-ray luminosity of approximately 37J 




(8) 



4'"" =^ 2.5 X 10^7 ergs-i 



(9) 
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• Fig. 4. Predicted polar cap heating luminosity, L+, for SCLF model vs. spin down 
age [37], Tsd, for different pulsar periods and an assumed polar cap surface tempera- 
ture of 3 X 10^ K. Thick curves are luminosities from curvature radiation (CR) pair 
heating and thin lines are luminosities from inverse Compton (ICS) pair heating. 
The solid circles show measured pulsar luminosities and upper limits of hot thermal 
components (see [45]). The dashed line is the fit of measured luminosity vs. age for 
the millisecond pulsars in the globular cluster 47 Tuc [291 lll| . 

Since nearly all millisecond pulsars produce only ICS pairs with incom- 
plete screening, the above expression overestimates the predicted L+ in these 
sources. 

4 Slot gap accelerator 

Due to the geometry of the field lines and the assumed boundary conditions of 
the accelerator, the altitude of the PFF varies with magnetic colatitude across 
the polar cap [5l[3S]. On field lines well inside the polar cap rim, is relatively 
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strong and the PFF is very near the neutron star surface. But at the polar cap 
rim, which is assumed to be a perfectly conducting boundary, E\\ vanishes. 
Near this boundary, the electric field is decreasing and a larger distance is 
required for the electrons to accelerate to the Lorentz factor needed to radiate 
photons energetic enough to produce pairs. The PFF thus curves upward as 
the boundary is approached, forming a narrow slot gap (see Figure [6|) near 
the last open field line [5]. Since E\\ is unscreened in the slot gap, particles 
continue to accelerate and radiate to high altitude along the last open field 
lines. The width of the slot gap is a function A = P ^2 of pulsar period 
and surface magnetic field f57], and can be expressed in magnetic colatitude 
as a fraction of the polar cap half-angle Zi^sG: where ^ = 9/dpc 

A. ^/4^, /I < 0.075 
"^^^^^(0.3, ^> 0.075 

The particles can achieve very high Lorentz factors which at altitudes of sev- 
eral stellar radii are limited by curvature radiation losses, to jsG — 3 — 4 x 10^ 
[58] . Since the slot gap is very narrow for young pulsars having short peri- 
ods and high fields, the corresponding solid angle of the gap emission fisG oc 
OpQrA^SG is quite small. So even though only a small fraction of the polar 
cap flux is accelerated in the slot gap, the radiated flux <1>SG = Lsg/^SG(P 
can be substantial. The total luminosity divided by solid angle from each pole 
is (from [57] ) 



^SG 
^SG 



e-y [0.123cos^a-|-0.516'ppsin^a] ergs"^sr"^ 

^ f 9 X 1034 L%\^PI'^ , B < 0. 1 Ser 

l2xl034Lf/^35P„«//, i3>0.1B„ 



(11) 



where Po.i = P/Q.l s, Lsd,35 = isd/lO^^ ergs^^ is the spin down luminosity, 
Bci = 4.4 X 10^3 G is the critical magnetic field strength and e-y is the efficiency 
of conversion of primary particle energy to high-energy emission. 

For A ^ 0.075, which corresponds to pulsars below the CR pair death 
line (see Fig. [2]), the slot gap disappears since the screening of is no longer 
effective. These "pair-starved" pulsars may have local screening of near 
the ICS PFF, but ^ at higher altitude so that acceleration can occur 
over nearly all of the open field volume [59]. The particle Lorentz factors, as 
in the slot gap, will be limited by curvature-radiation reaction to 7 ~ 10^. 



5 High-Energy Radiation 

5.1 Polar cap and slot gap cascades 

As discussed in Section 12.21 pair cascades above the polar cap may be ini- 
tiated by either ICS or CR. In the strong magnetic fields near the neutron 
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star, the Compton scattering cross section has a resonance at the cyclotron 
energy, ~ 12i?i2 keV, where its value is several orders of magnitude higher 
than the Thompson cross section. The thermal X-rays from the neutron star 
surface will be blueshifted into the resonance in the rest frame of the accel- 
erating primary electrons when kT"f{l — (3 cos 9) ~ es, where 9 is the angle 
between the photon propagation and the electron velocity (i.e. the magnetic 
field direction). Since the enhancement of the cross section at the resonance 
is so large, resonance scattering will dominate if the resonance condition is 
met [T51I75]. For typical pulsar fields and surface temperatures, electrons with 
7 ~ 10^ — 10^ will resonant scatter surface thermal X-rays [84l[35] to energies 
e ~ jB'mc^, where B' = B/Ba is the magnetic field strength in units of the 
critical field strength, Bcr = ^-^ x 10^'^ G. X rays scattered by electrons with 
7 ~ 10^ — 10^ will produce pairs to form a PFF. Resonant scattering is the 
dominant mode of ICS pair production for pulsars having the highest fields 
[30], whereas pulsars with lower fields, that include the millisecond pulsars, 
produce ICS pairs through non-resonant scattering. For 7 ~ 10^, the electrons 
will scatter non-resonantly in the Klein-Nishina limit, to energies e ~ 7. The 
ICS pair cascade multiplicities are Mi^^ ~ N+/Np ~ 10"'^ — 10 (pairs per 
primary electron) 40J, which is too low and produced at too low an altitude 
to completely shut off the acceleration. When the primary electrons reach 
7 ~ (0.5 — 1) X 10^, they produce CR photons that can produce pairs to form 
a PFF at height 0.02 — 0.1 stellar radii. The CR pair cascade multiplicities 
reach as high as M^^ ~ 10^ - 10"^ |40l [3]. 

The pairs are produced in Landau states whose maximum principal quan- 
tum number is rimax = 2e'(e' — 2)/B' |17| . where e' = esin^ is the photon en- 
ergy in the frame in which it propagates perpendicular to the local magnetic 
field. When the local field i?<0.1i?cr, pairs will be created above threshold 
in highly excited Landau states and the excitation level is quite sensitive to 
field strength, rimax cx 1/B'^. The pairs will decay through emission of syn- 
chrotron/cyclotron photons, many of which will produce more pairs in excited 
states. A pair cascade can be sustained in such a way through several gen- 
erations. The pairs can also produce ICS photons, through scattering surface 
thermal X-rays, which may produce pairs [82]. In high- field pulsars, where 
B ~ O.lBcr near the neutron star surface, the pair creation attenuation coeffi- 
cient is high enough for pair creation near threshold, so pairs are produced in 
very low Landau states. In this case, cascade pair multiplicities are lower since 
the number of cyclotron photons drop significantly • Bound pair production 
[68] also becomes important for B ^ 0.1i?cr, which will further lower the pair 
multiplicity. When B ^ i?cr, photon splitting dominates over pair production 
[33] . The screening of the electric field is also delayed by splitting and creation 
of bound pairs, effectively increasing the accelerating potential. The effect of 
bound pair creation on polar cap acceleration has been studied by Usov & 
Melrose [77]. 
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Radiated spectra from CR-initiated [Hj and ICS-initiated [75] cascades are 
very hard (roughly power laws with indices 1.5 - 2.0) [32] with sharp cutoffs 
due to magnetic pair production at energy [321 [30] 



where r is the emission radius and R is the neutron star radius. An approxi- 
mate form for the polar cap pair cascade spectrum is given by 



where a is the power-law index, Ci^ = Q.2{a / X){B' E? / p) and A — Xc/'^.ir is 
the electron Compton wavelength. The polar cap cascade radiation produces 
a hollow cone of emission around the magnetic pole, with opening angle, 
9ry ~ (3/2)(i7r/c)^/^, determined by the polar cap half-angle at the radius 
of emission r. The characteristics of emission from this type of polar cap 
model [19] has been successful in reproducing many features of 7-ray pulsars, 
including the wide double-peaked pulses observed from 7-ray pulsars like the 
Crab, Vela and Geminga, and the phase-revolved spectra. However, the polar 
cap opening angle is very small (a few degrees) unless the emission occurs 
more than a few stellar radii above the surface [Hj. Since the pair cascades 
over most of the polar cap occur within several stellar radii of the stellar 
surface, the broad profiles which require beam opening angles of the order 
of the magnetic inclination angle a, cannot be produced unless the pulsar is 
nearly aligned. Daugherty & Harding [TO] had to assume extended acceleration 
to 3 stellar radii and an artificial enhancement of primary particle flux near 
the polar cap rim in order to reproduce the Vela pulsar spectrum and pulse 
profile. 

Muslimov & Harding [S^ found that pair cascades on the inner edge of the 
slot gap occur at altitudes of 3-4 stellar radii and have higher multiplicities, 
M+ ~ 10^ — 10^, than the polar cap cascades. These are the characteristics 
required by [19] to model 7-ray pulsars, so a combination of polar cap cascades 
near the stellar surface and slot gap cascades near the polar cap rim at higher 
altitude may be successful in explaining non-thermal radiation from some 
fraction of sources having small magnetic inclination a and viewing angles C,. 
Figure [5] shows a plot of the total sky emission, phase vs. viewing angle 
C relative to the rotation axis, for a slot gap cascade at inclination angle 
a = 10°. In this model, the magnetic pole is at phase 0°, midway between the 
two peaks in the profile, which would be the predicted phase peak of polar 
cap thermal emission. 

5.2 Radiation from the high-altitude slot gap 

The electrons that accelerate in the slot gap and generate pair cascades at low 
altitude continue to accelerate. They will radiate curvature, inverse Compton 




(12) 




(13) 




Fig. 5. Phase plot, sample lightcurve, and a sketch of the accelerator location for 
the polar cap model, for a typical inclination angle a = 10° . The central zoom gives 
the gap extent relative to the star size. The dashed lines outline the null surface. 
The shading in the lightcurve and gap sketch is the same. The phase plot illustrates 
the change in lightcurve as seen by different observers and the aperture of the pulsed 
beams. From 128] . 



and synchrotron radiation at high altitudes. Initially, their Lorentz factors 
will be limited by curvature-radiation reaction, to 

/op ^2 \ 1/4 

7c«. = (^^) -3x10^ (14) 

and the peak energy of their CR spectrum will be Epeafc = "^^cIcrr/ P'^ ~ 
30GeV. As the electrons reach higher altitude, where the local magnetic field 
has dropped to B ^ 10^ — 10^ G, they may be able to resonantly absorb radio 
photons of energy eq^ghz that are at the cyclotron resonance in their rest 
frame . The condition for resonant absorption is 

^« = 3 X 10^- f ' (15) 

where Bs = B/10^ G and fiQ is the cosine of the angle between the radio 
photon and the electron momentum. The electrons are then excited to higher 
Landau states and radiate synchrotron emission. Blandford & Scharlemann 
[lOj investigated this process for pulsars, concluding that it was not an impor- 
tant source of radiation. However, they assumed that the electrons would be 
excited to only low Landau states and based their estimate on the cyclotron 
emission rate. Decades later, Lyubarsky & Petrova [51] re-examined resonant 
cyclotron absorption of radio emission in pulsar magnetospheres and discov- 
ered that the rate of absorption is much higher than the rate of re-emission 
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in low-lying Landau states, so the electrons will continue to be excited to 
high Landau states, and their pitch angles will increase until their momen- 
tum perpendicular to the magnetic field is relativistic. The emission is then 
synchrotron, not cyclotron, which proceeds at a much higher rate. Petrova 
[64] concluded that this process could be a significant source of radiation for 
young pulsars. It could also be significant for millisecond pulsars with high 
radio luminosities hke PSR J0218-I-4232 ([39] and Section 0|. 

In the slot gap or along open field line of pair-starved pulsars, electrons can 
experience both continuous acceleration and resonant cyclotron absorption. A 
steady-state may be reached between synchrotron radiation (SR) losses and 
absorption, with perpendicular momentum 

„SRR. 

^ — ~ 302B^'^Ey^ (16) 

which can indeed be relativistic. In this state, the Lorentz factor can remain 
locked at <IC Jcrr cyclotron resonance as the electron moves along 

the field with the field decreasing [39] . The synchrotron peak frequency de- 
creases with increasing altitude, spreading the emission over a range of ener- 
gies up to 100 MeV, and the curvature radiation rate drops dramatically. 

Non-resonant inverse Compton scattering of the radio emission by rela- 
tivistic electrons in the slot gap is also possible. The scattering would be 
completely in the Thompson regime, and the spectrum would extend to a 
maximum energy of Emax ~ 1crr^o,ghz ~ a few GeV. This process might 
make a significant contribution to the total spectrum of the slot gap if the 
radio emission region is located at relatively high altitude (see Section |6T]) . 



5.3 Relativity, geometry and caustics 

The geometry of emission at high altitude is strongly influenced by special rel- 
ativistic effects of aberration, time-of-flight and retardation of the magnetic 
field. Aberration and time-of-flight produce phase shifts of comparable mag- 
nitudes ~ —t/Rlc in radiation emitted at different altitudes r relative 
to the light cylinder radius Rlc- Morini [56' first noted that the combined 
phase shifts from aberration and time-of-flight, of photons radiated tangent 
to a magnetic dipole field from the polar cap to the light cylinder, nearly 
cancel those due to field line curvature on the trailing edge of the open-field 
region. Radiation along such trailing field lines bunches in phase, forming a 
sharp emission peak or caustic in the phase plot or profile (Figure (5]) . On the 
leading side, these phase shifts add up to spread photons emitted at various 
altitudes over a large range of phases. The effect is most pronounced for large 
a and emission between altitudes rem ~ 0.2 — 0.8 Rlc- Sweepback of the mag- 
netic dipole field due to retardation [20l [80] affects photon emission directions 
near the light cylinder, and also distorts the polar cap and open field volume 
[2l[22], even at the stellar surface. Most observer angles sweep across caustics 
from both magnetic poles, resulting in a double peaked pulse profile where the 
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Fig. 6. Same as Figure [5] for the slot gap model, for a typical inclination angle 
a = 45°. From [28]. 

peaks generally have phase separation less than 180° [24 . Furthermore, emis- 
sion occurs at all phases in the profile. Such profiles are very similar to those 
of the bright 7-ray pulsars. Crab and Vela. The predicted polarization char- 
acteristics of such a "two-pole caustic" model can also explain the observed 
optical polarization of the Crab pulsar [23] . 

Radiation from the slot gap has a geometry very similar to that of the 
"two-pole caustic" model and thus displays caustics in the intensity phase 
plots as well as Crab-like pulse profiles [SSj- The high-altitude slot gap thus 
may be a viable model for high-energy emission from young pulsars. However 
particle acceleration in the slot gap may not operate at all inclination angles. 
The low altitude E^\, shown in Eqn ([5]), is a function of both a and (p. For 
large a and cost/) < 1 (field lines curving away from the rotation axis), the 
second term in Eqn ([5]) can become large and negative at large r, causing the 
ii'll to reverse sign. This would cause a buildup of charge at that location, and 
would clearly be an unstable situation. The charge flow along those field lines 
might then be either time-dependent or non-existent. The solution for at 
high-altitude in the slot gap [58] may moderate this effect somewhat, but for 
fast rotators and a ~ 70° the reverses sign at low altitude. The resolution 
to this problem is still forthcoming, but it is possible that steady current flow 
in the slot gap does not occur for all geometries. 

5.4 Radiation from millisecond pulsars 

There are more than 200 radio pulsars now known with periods between 1 and 
30 milliseconds. Many (around 30) have been detected as X-ray point sources, 
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Fig. 7. Model spectrum for 2.3 ms PSR J0218+4232 (from 39 ) showing components 
from synchrotron radiation (SR) and curvature radiation (CR) from unscreened 
acceleration of a single primary electron, from the NS surface to the hght cyhnder, 
along a field line defined by magnetic colatitude ^ = 0.9 in units of PC half angle, 
at inclination a = 50°. Data points are from [48) . 



and seven of these have X-ray pulsations. Only the few millisecond pulsars 
(MSPs) that lie above or near the CR pair death line are expected to have slot 
gaps. These include B1821-24, B1957+20 and B1937+21, whose pulse profiles 
interestingly show caustic-like narrow, sharp peaks. The rest that lie below 
the CR death line are pair-starved and their SCLF is unscreened. Particles 
on all open field lines will therefore accelerate to high-altitude with Lorentz 
factors limited by CR losses, as in the slot gap (Eqn [13]). The peak energy of 
the CR spectrum for MSPs is similarly high, e^^^^ « XQQeN B^J ^ P~l' ^ kJI^, 
where Pms = P/lms and -Bg = G. But the very low magnetic fields of 

MSPs do not absorb the high energy part of the spectrum, so these sources 
may be visible up to energies of 50 GeV [TH [551 [35] . The CR will dominate for 
viewing angles near the magnetic poles because the emission comes primarily 
from low altitude, since the particle Lorentz factors decrease with increasing 
altitude. Additionally, as discussed in Section [5. 2[ synchrotron radiation from 
resonant absorption of radio emission may be important for MSPs having high 
radio luminosities. The synchrotron component will appear in the X-ray to 



Crab J1513-5908 Vela J1709-4429 J1952+3252 Geminga J1057-5226 



Radio 



Optical 



Soft 
X-Ray 



Hard 
X-ray 



y-ray > 
100 MeV 



y-ray > 
5 GeV 



Pulse Phase 

Fig. 8. Profiles of seven 7-ray pulsars in six energy bands, as in [73] . 

7-ray ( < 100 MeV) region of the spectrum [5^1, and will be visible at a larger 
range of observer angles since the emission comes from high altitudes. Figure 
[7] shows a model spectrum for the millisecond pulsar PSR J0218+4232 based 
on a (ID) cyclotron resonant absorption model described above (see Section 
lOl) . 



6 Pulsar Emission at Multiwavelengths 

Since rotation-powered pulsars shine over a broad spectrum from radio to 
high-energy 7-ray wavelengths, the multiwavelength spectra and profiles can 
give important clues to the acceleration and emission geometry. Although 
there are some clear patterns of spectral behavior with pulsar age [9], the 
observed multiwavelength lightcurves show a wide variety of characteristics 
that present many puzzles. As shown in Figure[8l the lightcurves for the known 
7-ray pulsars do not as a rule exhibit much phase correspondence between 
radio. X-ray and 7-ray peaks. The sole example is the Crab pulsar, where two 
peaks separated by ^ 140° are in phase from radio to 7-ray energy. For the 
other pulsars, the radio peak usually leads the one or two 7 peaks in phase and 
the soft X-ray peaks are broader and overlap the 7-ray peaks. Comparing X- 
ray and radio profiles of known rotation-powered X-ray pulsars, one sees that 
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there is more multiwavelength phase coherence for shorter periods {P < 50 
ms), such as for PSR B0540-69 (50 ms), PSR J1617-5055 (69 ms) and the 
miUisecond pulsars, PSR B1821-24, PSR B1937+21 and PSR J0218+4232. 
Such short-period and younger pulsars also tend to have pure power law X-ray 
spectra with no thermal components [9] , indicating that their X-ray emission 
is primarily non-thermal radiation from accelerated particles. One possible 
picture is that pulsars with shorter periods accelerate particles mostly in the 
outer magnetosphere, in slot gaps or outer gaps. The high-energy emission 
occurs at high altitudes that are large fractions of their light cylinder radii, 
i?/,^ = 5 X lO^P^^ cm, which are also smaller. Their emission geometry is 
therefore expected to be different from that of longer period pulsars, whose 
high-energy emission occurs closer to the neutron star surface and at altitudes 
that are very small fractions of Rlc- 

6.1 Radio emission geometry 

Because the mechanism responsible for the radio beams is not understood, and 
more importantly because the radiation is coherent, it has not been possible 
to describe this emission using a physical model. The emission has therefore 
been described using empirical models, developed over the years through de- 
tailed study of pulse morphology and polarization characteristics. The emis- 
sion is also highly polarized, and displays changes in polarization position 
angle across the profile that often matches the swing expected for a sweep 
across the open field lines near the magnetic poles in the Rotating Vector 
Model t65j. Empirical models (e.g. [66]) characterize pulsar radio emission as 
having a core beam centered on the magnetic axis and one or more hollow 
cone beams outside of the core. The average-pulse profile widths and compo- 
nent separations are measured to be decreasing functions of radio frequency 
([69]). This is consistent with a hollow cone beam centered on the magnetic 
pole, emitted at an altitude that decreases with increasing frequency (radius- 
to-frequency mapping). From such profile width measurements, and assuming 
that the edges of the pulse are near the last open field line, Kijak & Gil [46] 
find an emission radius of 



P<0.1 s the emission radius rradio ^ 0.05 — 0.1 Rlc- Such a picture is inde- 
pendently supported by polarization studies [A2] of pulsars younger than 75 
kyr. They concluded that the emission of these pulsars was from a single wide 
cone beam, that core emission was weak or absent, and that the height of the 
cone emission is between 1% and 10% of the light cylinder radius. 




(17) 



where P is the period derivative. This result predicts that for pulsars with 
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6.2 The global picture 

As discussed in Section [573l emission along the last open field lines at altitudes 
fern ^ 0.2 — 0.8 Rlc will form caustics, producing sharp peaks in the pulse 
profiles. According to Eqn (fT7|) . the radio cone emission altitude for pulsars 
with periods P < 30 — 50 ms, depending on their P, will fall in the range of 
caustic formation. If the radio conal emission is radiated over an extended 
range of altitude (i.e. a few tenths of Rlc) then radio caustic peaks would 
appear in phase with high-energy caustic peaks. Such a model would explain 
the multiwavelength phase coherence of the profiles of fast pulsars like the 
Crab and millisecond pulsars, PSR B1821-24 and PSR 61937+21. The conal 
emission of slower pulsars would fall below the altitude range of caustic forma- 
tion and the radio conal peak(s) or core emission would lead the high-energy 
caustic peaks, as seen in most of the multiwavelength profiles of Figure [S] In 
these cases, a single radio peak leading double 7-ray peaks is most likely to 
be the edge of a cone beam, since it is unlikely that a viewing angle crosses 
the narrow core beam as well as both caustics. For example, the observer 
viewing angle illustrated in Figure [H would cross the outer edge of pole 2 at a 
phase of 180°. Due to the radius-to-frequency mapping of the cone emission, 
this picture would predict that the measured width of the radio beam would 
be smaller at higher frequency since the emission is originating at lower al- 
titude, whereas the core emission width would not be expected to vary with 
frequency. 

Aside from the geometry of the emission at different wavelengths, that is 
directly relevant to observational characteristics of pulsars, there is the more 
fundamental question of how the microphysics of acceleration and pair cas- 
cades fit into the studies of the pulsar magnetosphere. The magnetosphere 
models [HI [71] make the assumption of ideal MHD (no parallel electric fields) 
in order to derive the global structure of the magnetic field and currents that 
are solutions to the so-called pulsar equation [54] that describes a spinning 
neutron star with a dipole field. A notable feature of these models is the 
formation of a neutral sheet in the equatorial plane of the spin axis which 
may provide a source of the return current that flows back to the neutron 
star along the last open fleld line. The derived currents in such models do 
not match the currents that have been assumed in either polar cap or outer 
gap accelerator models. Yet the microphysical acceleration models are needed 
to describe the physics of the pair creation and to check the consistency of 
the assumption of ideal MHD conditions over most of the magnetosphere. At 
present, it is not clear whether the global models can be adjusted to match the 
boundary conditions of the acceleration models [75] or the whether the accel- 
eration models can produce the required cascade multiplicity with compatible 
boundary conditions. In the end, self-consistency between global models and 
acceleration models may require time-dependence or spatial-dependence of 
the current flow. 



20 Alice K. Harding 



7 Open Questions 

Although much of the theoretical picture presented in this article depends on 
a set of undisputed fundamentals, such as, that rotation-powered pulsars are 
neutron stars with large-scale dipole magnetic fields that are directly measur- 
able from their period derivatives, to within uncertainties in the equation of 
state. But the picture also involves a number of underlying assumptions that 
are less well confirmed or agreed upon. The most important of these for par- 
ticle acceleration are the boundary conditions on the charge flow. Both polar 
cap and outer gap acceleration models assume that the boundary between the 
open and closed field line region is a perfect conductor. This implicitly assumes 
that there is some microphysical screening at this boundary by charges in the 
closed-field region, although this process has not been modeled or simulated. 
The boundary condition on the electric field and potential at the neutron star 
depends on the binding energy of charge in the surface layers. Calculations 
of the cohesive force for a solid lattice in magnetic fields up to about 10^'^ G 
indicate that charges will not be bound in the surface at temperatures above 
~ 0.5 MK, so that most pulsars have surfaces hot enough to allow SCLF 
accelerators to operate. However, a number of pulsars have magnetic fields 
above 10^'^ G, where surface conditions have not been adequately explored. It 
is possible that the neutron star surface may be in a liquid rather than solid 
state in very high magnetic fields ([43l [34]) or that charges may be bound 
by much stronger forces. In this case, high-field pulsars could have vacuum 
instead of SCLF accelerators. 

Models often also assume that the neutron star is a dipole field in calcu- 
lating the accelerating electric potential and field. This assumption could be 
faulty at low altitude if there exist higher multipole fields near the neutron 
star surface, and is certainly invalid at very high altitude where there are 
distortions of the field due to retardation [22] and current flow ([60]). The E'n 
in the high altitude slot gap [5S] is only approximate and needs to be calcu- 
lated using field line distortions. It is possible that the distorted field lines, 
having less curvature than a dipole, will have the effect of removing the sign 
reversals present in the approximate solutions. There are also expected to be 
cross field particle motions at high altitude that will modify solutions of the 
electric field. 

The possible connection of pair cascades to radio emission morphology 
is an intriguing question that needs further investigation. For example, are 
the pairs from the near-surface polar cap and high-altitude slot gap cascades 
responsible for the core and cone beams? Do pairs produced in the outer gaps 
produce any radio emission? Whether polar caps and slot gaps can exist on 
the same field lines or in the same pulsar magnetosphere is also not known. 

Although some of these questions require more theoretical work, it is 
hoped that many answers will come from observations with future detectors. 
The Large Area Telescope (LAT) on the Gamma-ray Large Area Space Tele- 
scope (GLAST), due to launch in 2007, will have unprecedented sensitivity 
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and energy resolution for gamma-rays in the range of 30 MeV to 300 GeV. 
GLAST is therefore expected to provide major advances in the understand- 
ing of high-energy emission from rotation-powered pulsars [741 131j , including 
an increase in the number of detected radio-loud and radio-quiet gamma- 
ray pulsars, and millisecond pulsars, giving much better statistics for eluci- 
dating population characteristics, measurement of the high-energy spectrum 
and the shape of spectral cutoffs and determining pulse profiles for a variety 
of pulsars of different ages. Further, measurement of phase-resolved spectra 
and energy dependent pulse profiles of the brighter pulsars should allow de- 
tailed tests of magnetospheric particle acceleration and radiation mechanisms. 
Third-generation ground-based Air Cherenkov detectors, such as H.E.S.S. [61] 
in Namibia, MAGIC [50] in the Canary Islands and VERITAS [41] in the US, 
have begun operation. They are sensitive to 7-rays in the range 50 GeV to 
50 TeV and are putting important constraints on pulsar spectra and emission 
mechanisms [Tj. Further into the future are planned X-ray telescopes, such as 
Constellation X, XEUS, and X-ray polarimeters, such as AXP, POGO and 
ACT. Polarimeters [78] in particular will be extremely important in verify- 
ing model predictions for different pulsar emission geometries through phase- 
resolved measurements of position angle and polarization percent. 
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